Exact Ground State Properties of Disordered Ising-Systems 
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Exact ground states are calculated with an integer opti- 
mization algorithm for two and three dimensional site-diluted 
Ising antiferromagnets in a field (DAFF) and random field 
Ising ferromagnets (RFIM), the latter with Gaussian- and 
bimodal-distributed random fields. We investigate the struc- 
ture and the size-distribution of the domains of the ground 
state and compare it to earlier results from Monte Carlo sim- 
ulations for finite temperature. 

Although DAFF and RFIM are thought to be in the same 
universality class we found essential differences between these 
systems as far as the domain properties are concerned. For 
the DAFF the ground states consist of fractal domains with a 
broad size distribution that can be described by a power law 
with exponential cut-off. For the RFIM the limiting case of 
the size distribution and structure of the domains for strong 
random fields is the size distribution and structure of the clus- 
ters of the percolation problem with a field dependent lower 
cut-off. Consequently, the domains are fractal and in three 
dimensions nearly all spins belong to two large infinite do- 
mains of up- and down spins - the system is in a two-domain 
state. The fractal dimensions for the DAFF and the RFIM 
agree. The DAFF ground state properties agree with results 
from MC simulation in the whole whereas there are essential 
differences between our exact ground states calculations and 
earlier MC simulations for the RFIM. 

PACS: 75.50.Lk, 64.60.Cn, 75.40Mg 

Keywords: Ising-Models, Random Magnets, Numerical Meth- 
ods 

I. INTRODUCTION 

Many aspects of the influence of random field disorder 
on a system of interacting spins are still not well under- 
stood (for a review see [0 ) . Since it has been argued that 
the diluted Ising antiferromagnet in a homogenous field 
(DAFF) and the random field Ising ferromagnct (RFIM) 
are in the same universality class 
vestigations focus on the DAFF | 
realizations of a system with random field disorder while 
theorists usually focus on the RFIM Q . 

However, there have been investigations on the do- 
main structures of the RFIM [§ and the DAFF 
both based on MC simulations which suggest that there 
might be essential differences between DAFF and RFIM, 
at least as far as the domain structure is concerned. In 
the limit of strong disorder for the RFIM the domains 
are fractal in three and compact in two dimensions while 
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the domains of the DAFF have a fractal structure both 
in two and three dimensions. Additionally, due to the 
fractality of the domains essential deviations from Imry- 
Ma-type ^ behavior which is thought to be valid for 
small disorder have been found. 

In order to prove how far the results from MC sim- 
ulations which are out of equilibrium for these systems 
with frozen dynamics | [lO| can be transferred to the equi- 
librium T — case we performed exact ground state 
calculations on DAFF and RFIM with an integer opti- 
mization algorithm. With this method we get exact in- 
formation on equilibrium properties at zero temperature. 
We investigate the distribution of domain sizes and the 
structure of the domains and we compare our findings 
for the DAFF and the RFIM. Hopefully, the results may 
help to understand the critical behavior of these systems 
as well as dynamical aspects [|ll|. 

The Hamiltonian of the DAFF in units of the nearest 
neighbor coupling constant J is 



H 



E 

<ij> 



(1) 



with the uniform field i? > on all sites of the quadratic 
L X L respectively cubic L x L x L lattice. Here ai = ±1 
denotes Ising-spins and a fraction p of the sites is oc- 
cupied with a spin (quenched disorder: = 0,1). The 
phase diagram of the 2d-DAFF consists of an antiferro- 
magnetic low temperature phase for magnetic field B = 
and a disordered phase for all finite values of -B. In 
three dimensions there exists a long-range ordered phase 
Jl2| , p^ for magnetic fields B smaller than the critical 
fields Be which is Be « 1.4 for a dilution of p = 0.5 
0. For higher fields and low temperatures the DAFF 
develops a frozen disordered domain state in two and 
three dimensions. This domain state has many of the 
characteristics of a spin glass, as for instance a remanent 
magnetisation and an irreversibility line scaling like the 
deAlmeida-Thouless line . 

The Hamiltonian of the RFIM in units of the nearest 
neighbor coupling constant J is 



H 



E 

<ij> 



E^> 



(2) 



Here, all sites are occupied and the random fields Bi are 
taken from a Gaussian (Gaussian-RFIM) respectively bi- 
modal (±A-RFIM) probability distribution correspond- 
ing to 
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and 

P(i?,)-i(<5(B,-A)+<5(i?, + A)), (4) 
respectively. 

As in the case of the DAFF the phase diagram of the 2d 
system consists of a long-range ordered low temperature 
phase for zero random fields and a disordered phase for all 
finite values of the random field ||l^ . In three dimensions 
for weak random fields and low temperatures there exists 
a long-range ordered phase |]l3|. The critical value for 
the random field at which long-range order disappears 
was found to be Ac = 2.35 at zero temperature [ p^ . 
Even for the RFIM a possible glassy transition has been 
discussed [ p^ . 

We calculate exact ground states using an optimization 
algorithm well known in graph-theory. The Ising-system 
is mapped on an equivalent transport network, and the 
maximum flow is calculated using the Ford-Fulkerson al- 
gorithm [p6|-p^ . In order to investigate the distribution 
and the structure of domains we perform a cluster anal- 
ysis with a suitable adjusted Hoshen-Kopelman type al- 
gorithm |]l^. This algorithm pieces the system into do- 
mains, where a domain is defined as a group of spins 
which are connected and antiferromagnetically (DAFF) 
respectively ferromagnetically (RFIM) ordered. 

The investigations for the DAFF are carried out for 
systems of size 400 x 400 with an average over 50 different 
dilution configurations with p = 0.7 for each value of B. 
For some structural aspects even ground states of size up 
to 700 X 700 were calculated. Systems of size 50 x 50 x 50 
averaged over 40 dilution configurations with p — 0.5 are 
used for = 3. The calculation of ground states for the 
RFIM is more time consuming. Here we considered 300 x 
300-systems with 50 different random field configurations 
in 2d and 40 systems of size 30 x 30 x 30 in 3d. 

II. DISTRIBUTION OF DOMAIN-SIZES 



domain sizes. This relation indicates that except from 
the isolated spin-clusters only large domains can be real- 
ized for very weak fields. Due to this finite size effect we 
consider here fields i? > 0.5 only and hence, all domain 
sizes Vd > 3 exist. 

As has been shown earlier for the two dimensional 
case the distribution of domains in a DAFF is well de- 
scribed by a power law with an exponential cut-off, 

NDiVo) ~ V^^ eM-Vo/Vo). (5) 

where Vq depends strongly on the field. Because of the 
inaccuracy in the determination of Vq no significant p- 
dependence can be observed for this quantity. We in- 
creased the numerical effort in order to investigate a 
possible concentration and field dependence of S. We 
find S weakly decreasing for increasing field but no p- 
dependence. The results are summarized in Tab. ^ 

Fig. shows the size-distribution of the domain state 
of the three dimensional DAFF for different fields. For 
the largest value of the field a reasonable fit to Eq. |^ can 
be carried out yielding S — 1.8 ± 0.3 and Vq = 50 ± 20. 
However, a detailed investigation of the field or concen- 
tration dependence of Vq and S is not possible with sat- 
isfactory precision for the following reasons: As has been 
shown earlier |^ for B « Be the domain state consists 
mainly of two large interpenetrating domains (for B = A 
these domains contain approximately 80% of the spins) . 
In Fig. the two isolated peaks in the size distribution 
corresponding to the two percolating domains are not 
shown but due to the existence of these two infinite do- 
mains there are only few data for domains of finite size so 
that the statistics for these domains is rather bad. In the 
limit of very high fields the domains become very small. 
For S > 6 all spins are polarized by the field, i. e. all 
domains have size 1. Hence, there is only a small region 
of values of the magnetic field where a broad distribu- 
tion of domain sizes exist. As we will show in the next 
subsection the existence of two infinite domains has even 
more dramatic consequences for the RFIM. Here, there 
are no values of random fields where relevant domains of 
finite size can exist. 



A. DAFF 

In two dimensions for B > no long-range order exists 
and the ground state of the system consists of antiferro- 
magnetic domains of finite size. Within a ground state for 
each domain the domain- wall energy E^, i. e. the number 
of broken bonds, must be smaller than the field-energy 
BrriD, where m/^ is the absolute value of the magnetisa- 
tion of the domain which is not zero in these antiferro- 
magnetic domains due to the dilution. Therefore, only 
domains with mu > Eii,/\B\ exist in the system. The 
minimum wall energy of a domain is one broken bond, 
= 1. Additionally, it is always Vd > mo, where Vd 
is the size of the domain (number of spins), so that with 
Vd > we get a rigorous lower cut-off for possible 



B. RFIM 

In the limit of strong random fields A > 2d (±A- 
RFIM) respectively A ^ oo (Gauss-RFIM) aU spins fol- 
low the direction of the random field. Hence, the fer- 
romagnetic domains are the clusters of sites with uni- 
form random field signs. In this limit the domains of the 
RFIM correspond to clusters of the percolation problem 
pO| with a concentration of 50%. 

The distribution of domain-sizes for different A is 
shown in Fig. ^ for the two dimensional ±A-RFIM. For 
comparison the size-distribution for the percolation prob- 
lem with a concentration of 50% is also shown. With in- 
creasing A, the domain-size distribution approaches the 
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cluster-size distribution of the percolation problem which 
follows Eq. I with S = 1.55 ± 0.05 and Vq = 120 ± 20. 

For lower random fields there is a striking difference 
between the Gaussian-RFIM on the one hand and the 
±A-RFIM on the other hand: for the latter system there 
is a lower cut-off for the possible domain-sizes depending 
on the strength of the random field A. This cut-off can 
be seen in Fig. I and it can be determined analytically: 
Similar to the arguments given for the DAFF for each 
domain the relation < Eb — 'Y^=i^i^i niust hold. 
The positive domain-wall energy has to be compensated 
by the negative field energy. In the most favorable case 
all random fields acting on the spins of a domain have the 
same sign, Eb = VdA so that for the smallest possible 
domains we get Vd > Eyj/ A.. With this relation as start- 
ing point a lower cut-off can be determined as the sizes of 
those domains which have a minimum surface. For large 
domain sizes the minimum surface of a domain is the 

1/2 

surface of a circle (2d), Ey^ ~ Vjj , respectively sphere 
(3d), E^ - Vl'^. For smaller sizes, the domain with the 
minimum surface Vmm can be determined numerically. 

Fig. H shows this dependence of Vmin from A for the 
±A-RFIM in two dimensions. In addition, the minimum 
domain-volumes of the ground states of the RFIM for 
different A are depicted. For A < 1 there is a strong 
deviation of the observed domains of minimum size from 
the theoretical curve. Obviously, this is due to the fact 
that for larger domains it is less probable that all ran- 
dom fields within a domain with minimized surface have 
the same orientation. Hence, either the field energy Eb 
is overestimated or the shapes of the domains are more 
complicated than those with minimum surface. The first 
possibility would correspond to a crossover to a behav- 
ior following the Imry-Ma-argument [H, where Eb is 
thought to scale with the root of the domain size due 
to simple statistical fluctuations of random fields within 
a domain. However, in the range of fields where this 
crossover may occur the minimum domain size is larger 
than the system sizes that we can investigate with nu- 
merical methods. 

Fig. ^ shows the volumes of the two largest domains 
in the ground state of the ±A- and the Gaussian-RFIM. 
For lower fields A < 0.8 the order parameter is finite. 
Obviously, this is the region where finite size effects lead 
to wrong results for the system sizes we investigated since 
in this regime minimum domain sizes are larger than the 
system. Therefore, in the following we consider only re- 
sults for A > 0.8. 

An additional, striking feature of Fig. ^ is the jump 
in the data for the ±A-RFIM at A = 2. Considering 
Fig. ^ (left) one realizes A = 2 as a special value: The 
random fields are marked with -I- and — signs and the ar- 
rows represent the spin-orientations (t-spins are oriented 
parallel to a positive A). The energy of the center spin 
with cr = ±1 is 

E„^2a- a A. (6) 



In the ground state it is tr = -1-1 for A > 2 and the bold 
drawn domain wall is realized. A < 2 favors a = —1 and 
the center spin belongs to the down-domain. Generally, 
for A > 2 the domain wall arranges in such a way that 
the sign of the random fields changes at each position of 
a domain wall, i. e. the domain walls run along the sur- 
faces of the random field clusters (clusters of sites with 
uniform random field sign). Fig. ^ (right) demonstrates 
that nevertheless even for A > 2 domains larger than 
these random field clusters can exist. A cluster with a 
magnetisation that is oriented antiparallel to the random 
fields is enclosed such that its surface lies within the do- 
main. Hence, for A > 4 each domain is a random field 
cluster while for 2 < A < 4 each domain wall is a part of 
the surfaces of the random field clusters. Consequently, 
the domains in this range of A are larger. In three di- 
mensions each spin has six next neighbors. Here, the 
corresponding effect occurs at A = 4. 

Both, for d — 2 and d — 3 the RFIM approaches the 
percolation problem with a concentration of 50% for large 
A, but while for d = 2 in this case there exists no perco- 
lating cluster of random fields in three dimensions 50% 
is above the percolation threshold. This means that as 
mentioned above for the case of the DAFF in three di- 
mensions one will always find two infinite interpenetrat- 
ing domains within the disordered phase whose antipar- 
allel oriented magnetisations cause the order-parameter 
to vanish. The normalized volumes of the largest and 
second-largest domain for the d = 3 ±A-RFIM are shown 
in Fig. ||. Adding the sizes of these two largest domains 
we found that more than 97% of the spins of the system 
belong to these two domains in the whole range of ran- 
dom fields. In contrast to d = 2 there is no broad range 
of domain-sizes for d = 3 and in contrast to the DAFF 
there is no range of fields where a relevant number of do- 
mains of finite size exist. For this reason it is senseless to 
analyze the distribution of the domains in more detail. 

III. DOMAIN STRUCTURES 
A. DAFF 

In a DAFF domains are composed in such a way that as 
many spins as possible are oriented parallel to the applied 
field with a minimum number of broken bonds. Hence, 
domain walls favorably run along vacancies to minimize 
the wall energy. These two competing requirements for 
minimizing the total energy lead to highly non-trivial 
structures. Fig. |^ shows the antiferromagnetic domains 
in the ground state. Here, the configuration of occupied 
sites and vacancies is kept fixed, while the external field 
B is tuned. In order to investigate the structural proper- 
ties of the domains we consider the following quantities 
and corresponding scaling relations: 

• volume Vd- number of spins 
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surface Sd- sum of all unsatisfied bonds (broken or 
"open" due to vacancies) with Sd ~ Vj 



• radius of gyration Ru, i. e. root of mean squared 
distance of all spins of a domain with Vd ^ 



• absolute value of the magnetisation of a domain 



niD with wiD ^ V'j 



• wall energy Ew, 
Ew 



D 



sum of all broken bonds with 



For large domains the scaling relations above hold and 
the exponents are determined by analyzing all domains 
in the ground states and taking an average of the quan- 
tities above for all domains of equal size Vd- In we 
showed an example for the two dimensional case. In the 
meantime we investigated the three dimensional system 
also. Our resulting exponents are summarized in Tab. || 
and compared to the values resulting from MC Simula- 
tions 

The domains are fractal with fractal dimension = 
1.65 for d = 2 and D„ = 2.14 for d 3. This frac- 
tality is also reflected in the proportionality of surface 
and volume {Dyf « 1). A great part of the surface is 
inside the domains. Slight deviations of exponents from 
exact ground states from those resulting from MC sim- 
ulation are due to a systematic inaccuracy which is not 
reflected by the error-bars representing the fluctuations 
of the data only. Probably, a small, systematic error is 
due to the fact that the data are slightly curved in a log- 
log plot so that the linear range for large domains can 
hardly be determined (see also the corresponding Fig. |l^ 
for the RFIM). From comparison with the corresponding 
Figures [QjH it is found that in spite of similar exponents 
the values md{VD) are larger for exact ground states - 
the domains of the ground states are obviously better 
optimized. 

In good agreement with MC simulations is « 1 
both in two and three dimensions. Altogether the frac- 
tal structures for strong disorder lead to deviations from 
Imry-Ma-type arguments |^ based on compact domains. 
Here, a statistically distribution of vacancies in a domain 
would yield w 1/2 which can clearly be ruled out by 
our data for the range of fields that we investigated. A 
field dependence of any of the exponents could not be 
found. 



B. RFIM 

As we showed in the previous section there is no broad 
distribution of domains for the three dimensional RFIM - 
domains of finite size are irrelevant for the ground state. 
Therefore, in this subsection we restrict ourselves to the 
case d = 2. 

Fig. 1^ shows the ferromagnetic domain configuration 
of a ±A-RFIM in the ground state. The configuration of 
signs of the random fields is kept fixed, while the value 



A of the field is tuned. As for the DAFF we study the 
structural quantities volume, surface and radius of the 
domains. Since there is no dilution the surface and the 
domain wall energy are identical. Also, the magnetisa- 
tion and the volume of the domains are identical. We 
study the following quantities: 

• local random field fluctuations h d , which is the ab- 
solute value of the sum of the random fields of a 
domain 

• random field fluctuations of all sites of a domain 
along the domain wall inside the domain 

• random field fluctuations of all sites of a domain 
along the domain wall outside the domain hout 

Fig. 1^ shows bin/So and bout/So each multiplied with 
the sign of the magnetisation of the considered domain 
for the ±A-RFIM. As we demonstrated with Fig. || for 
A > 2 the domain walls run exactly along the boundaries 
of random field clusters. Hence, in this range of fields it 
is \bin/SD\ = \bout/SD\ = 1- Even for smaller fields the 
domain walls are presumably located along these random 
field surfaces and for large domains these fluctuations are 
proportional to the surface of the domains. 

Corresponding to the DAFF the exponents D^, D^f 
and 9 with bo ^ are determined from appropriate 



fits (Fig. |1C|). In Tab. Ill the exponents are summarized. 
The exponents resulting from MC are from 

The investigations are carried out for relative strong 
random fields, A = 3.01 (Gaussian-RFIM) and A = 4.01 
(± A- RFIM), because then - as explained in the previous 
chapter - a broad range of domain sizes exists. It is found 
that the exponents do not depend on the kind of random 
field distributions within the given precision. 

In agreement with the findings for the DAFF the do- 
mains are fractal with nearly the same exponents Dy w 
1.6 and Dyf « 1. Also, it is 6* « 1 indicating that the 
random fields are distributed non-statistically within a 
domain. These results are in contradiction to the earlier 
investigations using MC methods where results were 
found corresponding to more compact domains. We in- 
terpret this discrepancy as being due to the fact that 
the MC simulations were carried out with A = 1 and 
T/ks = 1. For this value of A only very large domains 
can exist in the ground state. Therefore, the broad range 
of analyzed domains in results obviously from thermal 
fiuctuations. For larger fields and lower temperature we 
would expect agreement of these two methods. 

Note that all exponents found for the fractal behavior 
of the domains of the RFIM as well as the DAFF are very 
close to those for the lattice animals of the percolation 
problem |2^,|2^ . This is plausible since as mentioned be- 
fore the domains of the RFIM in the limit of high random 
fields have the same distribution and structure like the 
clusters from the percolation problem. As we discussed in 
the previous subsection systematic inaccuracies are not 
refiected by the error-bars and systematic error can be 
due to a slight curvature of the data. 
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Our results deviate once more from the assumption 
of Imry-Ma-type random field fluctuations for com- 
pact domain structures with 9 = 1/2 and bin/out ^ L^^^- 
However, no A-dependencc for our exponents can be de- 
termined as indication for a crossover to Imry-Ma behav- 
ior within the range of fields that we investigated. 

IV. CONCLUSIONS 

We investigated the distribution and structure of the 
domains of the ground state of DAFF and RFIM in two 
and three dimensions. Since except for the Gaussian 
RFIM there is always a lower cut-off for the minimum 
domain size one is restricted to rather large disorder due 
to finite size effects, i. e. high dilution and fields of the 
order of the spin-spin coupling constant. 

The domain-size distribution for the DAFF can be well 
described by a power law with a field dependent expo- 
nential cut-off for a broad range of the applied fields in 
2d and for strong fields {B > 3) in 3d. For the RFIM 
with increasing random field strength there is a continu- 
ous transition to the cluster-size distribution of the per- 
colation problem which also follows a power law with an 
exponential cut-off for both Gaussian- and bimodal dis- 
tribution of the random fields. For the latter a non triv- 
ial minimum domain-size is found, that can be estimated 
based on energetic considerations. 

Investigating structural properties we have found that 
the domains are fractal for the DAFF (d = 2, 3) as well 
as for the RFIM (d = 2) with the same fractal exponents. 
Surprisingly, the domain state of the RFIM in 3d is found 
to consist mainly (> 97%) of two infinite interpenetrating 
domains of opposite phase in the whole range of random 
fields for which the long range order is broken. In this 
sense, the phase diagram of the RFIM at T = consists 
only of two regions, a one-domain state (long-range order 
below the critical field) and this two-domain state. 

While for the RFIM the results for the structural scal- 
ing exponents differ essentially from earlier findings at 
finite temperature carried out using MC simulations, our 
results are in good agreement with Monte Carlo simula- 
tions for the DAFF both in two and three dimensions. 

The magnetisation of the domains and the distribu- 
tion of random fields within the domains, respectively, 
strongly deviates from the assumptions of the Imry-Ma 
argument which however can still thought to be correct 
in the limit of small disorder. A more detailed investiga- 
tion of the structure of domains for weak random fields 
is desirable. Unfortunately, for this the numerics exceeds 
current capacities. 
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D 

_D 





0.6 


1.76 ±0.06 


0.9 


1.76 ±0.06 


1.2 


1.68 ±0.04 


1.4 


1.61 ±0.06 


1.5 


1.58 ± 0.05 


2.2 


1.44 ±0.03 


2.5 


1.42 ±0.06 


2.8 


1.37 ±0.04 


3.5 


1.32 ±0.05 



TABLE I. Field dependence for the exponent S of the do- 
main-size distribution for the 2d DAFF 





MC Simulation 


Exact Ground States 


d 


2 


3 


2 


3 


e 


1.00 ±0.01 


0.996 ±0.001 


0.97 ±0.01 


0.99 ± 0.01 


a 


0.98 ±0.01 


0.995 ±0.001 


0.96 ± 0.02 


0.98 ± 0.03 




1.001 ±0.002 


1.000 ± 0.001 


1.00 ±0.01 


1.00 ±0.01 




1.56 ±0.03 


2.0 ±0.1 


1.64 ±0.06 


2.14 ±0.08 



TABLE II. Structural exponents for the DAFF 





MC Simulation 


Exact Ground States 




±A 


Gaussian 


±A 


Gaussian 








1.65 ±0.08 


1.67 ±0.08 




0.59 ±0.04 




0.98 ±0.05 


0.97 ±0.05 


e 


0.66 ± 0.04 




1.00 ±0.01 


1.00 ±0.01 



TABLE III. Structural exponents for the 2d RFIM; 
A — 4.01 for ±A- and A = 3.01 for Gaussian-distributed 
random fields 



FIG. 1. Nd versus Vd for the DAFF in three dimensions. 
For B = 4.01 a fit to Eq. | is also shown as solid line. 

FIG. 2. Domain-size distributions for the two dimensional 
±A-RFIM (Impulses) and size-distribution of the random 
field clusters (open symbols) 

FIG. 3. Minimum domain size for given values A of the 
field of a ±A-RFIM in two dimensions 



FIG. 4. Largest (open symbols) and second largest (filled 
symbols) domain size divided by the system size Lx L for the 
±A-RFIM in two dimensions 

FIG. 5. ±A-RFIM; left: cut of a domain wall; right: a 
domain for A > 2 consisting of multiple random field clusters 



FIG. 6. Largest and second-largest domain size corre- 
sponding to Fig. ^ for the RFIM in three dimensions. 

FIG. 7. Ground states of a 200 x 200 DAFF with a fixed 
dilution configuration for increasing fields B = 0.8, 1.4, 2.2 
(from above). The two antiferromagnetic phases are repre- 
sented black and white the vacancies are also black. 

FIG. 8. Ground states of a 200 x 200 ±A-RFIM for a fixed 
configuration of random fields and A — 1.0, 1.1, 1.5 (from 
above). The two ferromagnetic phases are represented black 
and white. 

FIG. 9. Random field fiuctuations bout and bi„ of the 
domain surfaces each multiplied with the sign of the do- 
main-magnetisation and divided by the surface Sd for the 
two dimensional ±A-RFIM. 

FIG. 10. Figures for the determination of the structural 
exponents 6, and D^f for the two dimensional ±A-RFIM. 
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